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ABSTRACT 

, MESS (Mass-loss of Evolved StarS) is a Guaranteed Time Key Program that uses the PACS and SPIRE instruments on board the 

. Herschel Space Observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and 

^S) ' post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. 

In total, of order 150 objects are observed in imaging and about 50 objects in spectroscopy. 

This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated 
using results obtained during HerscheVs. science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB 
stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type 
stars. 



Key words, circumstellar matter - infrared: stars ~ stars: AGB and post-AGB - stars: mass loss - Supernova remnants - planetary 
nebulae 

1. Introduction phase, while for massive stars (initial mass >15 M©) the mass 

, . , J . ^. „ . , ^ . loss takes place in a fast (hundreds to a few thousand km s^h 
Mass-loss IS the dommatms factor m the post-mam sequence , j - i- j i 

, ^. „ , ° „ , A ■ ^ J- wmd driven by radiation pressure on lines at a moderate rate of 

evolution of almost all stars. For low- and intermediate mass ^ j^q_6 j^-^ /p^Jg 2008) 

stars (initial mass <8 M©) this takes place mainly on the ^ yr ( u s e a . 

thermally-pulsing AGB (asymptotic giant branch) in a slow (typ- Although mass loss is such an important process and has 

ically 5-25 km s-^) dust driven wind with large mass loss rates been studied since the late 1960's with the advent of infrared 

(up to 10-4 Mo yr-\ see the contributions in the book edited by astronomy, many basic questions remain unanswered even after 

Habing & Olofsson 2003), which is also the driving mechanism important missions such as IRAS (Neugebauer et al. 1984), ISO 

for the slighdy more massive stars in the Red Supergiant (RSG) (Kessler et al. 1996), Spitzer (Werner et al. 2004) and AKARI 

(Murakami et al. 2007): what is the time evolution of the mass- 

* Herschel is an ES A space observatory with science instruments pro- loss rate, what is the geometi'y of the mass-loss process and how 

vided by European-led Principal Investigator consortia and with impor- does this influence the shaping of the nebulae seen around the 

tant participation from NASA. central stars of Planetary Nebulae (PNe) and Luminous Blue 
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Variables (LBVs), can we understand the interaction of these 
winds with the interstellar medium (ISM) as initially seen by 
IRAS (e.g. Stencel et al. 1988) and confirmed by AKARI (Ueta 
et al. 2006, 2008) and Spitzer (Wareing et al. 2006), what kind 
of dust species are formed at exactly what location in the wind, 
what are the physical and chemical processes involved in driv- 
ing the mass-loss itself and how do they depend on the chemical 
composition of the photospheres? With its improved spatial res- 
olution compared to ISO and Spitzer, larger field-of-view, better 
sensitivity, the extension to longer and unexplored wavelength 
regions, and medium resolution spectrometers, the combination 
of the Photodetector Array Camera and Spectrometer (PACS, 
Poglitsch et al. 2010) and the Spectral and Photometric Imaging 
Receiver (SPIRE, Griffin et al. 2010) observations on board the 
Herschel Space Observatory (Pilbratt et al. 2010) have the po- 
tential to lead to a significant improvement in our understanding 
of the mass-loss phenomenon. This is not only important for a 
more complete understanding of these evolutionary phases per 
se, but has potentially important implications for our understand- 
ing of the life cycle of dust and gas in the universe. 

Dust is not only present and directly observable in our 
Galaxy and nearby systems like the Magellanic Clouds, but is 
already abundantly present at very early times in the universe, 
e.g. in damped Lyman-alpha systems (Pettini et al. 1994), sub- 
miUimetre selected galaxies (Small et al. 1997) and high-redshift 
quasars (e.g. Omont et al. 2001, Isaak et al. 2002). The inferred 
far-IR (FIR) luminosities of samples of 5 < z < 6.4 quasars are 
consistent with thermal emission from warm dust (T < 100 K), 
with dust masses in excess of 10^ solar masses (Bertoldi et al. 
2003, Leipski et al. 2010). 

It has been typically believed that this dust must have been 
produced by core-collapse (CC) SuperNovae (SNe), as AGB 
stellar hfetimes (10*^ to 10^ yr) are comparable to the age of 
universe at redshift > 6 (Morgan & Edmunds 2003, Dwek, 
Galliano & Jones 2007). The observed mid-IR emission for a 
limited number of extra-galactic SNe implies dust masses which 
are generally smaller than 10"^ M© (e.g. Sugerman et al. 2006, 
Meikle et al. 2007, Blair et al. 2007, Rho et al 2008, Wesson et 
al. 2010a), corresponding to condensation efficiencies which are 
at least two orders of magnitude smaller than theoretical mod- 
els predict (Todini & Ferrara 2001, Bianchi & Schneider 2007). 
FIR and sub-mm observations of dust within supernova rem- 
nants (SNR) estimate masses ranging from 0.1-1 Mg (Duime 
et al. 2003, 2009, Morgan et al. 2003, Gomez et al. 2009), yet 
there are a number of difficulties with the interpretation of these 
results. It is obvious that there is now indeed clear observational 
evidence for dust formation in CCSNe, but the quantity of dust 
formed within the ejecta is still a subject of debate. Valiante et 
al. (2009) recently showed that AGB stars could potentially ri- 
val or surpass SNe as the main producer of dust at characteristic 
timescales of between 150 and 500 Myr, although the model re- 
quires rather extreme star formation histories, a top-heavy initial 
mass function and efficient condensation of dust grains in stellar 
atmospheres. The dust production of SNe, either from the pro- 
genitors (LBV, RSG, Wolf-Rayet (WR) stars) or direcdy in the 
ejecta, versus that of AGB stars is therefore of utmost impor- 
tance and one of the science themes that will be addressed in the 
MESS Herschel key program described in this paper. 

Most of the astronomical solid state features are found in 
the near-IR (NIR) and mid-IR (MIR) ranges. The ISO SWS 
and LWS spectrometers revolutionised our knowledge of dust 
and ice around stars. In the LWS range, partly overlapping with 
Herschel PACS, most of ISOs spectroscopic dust observations 



suffered from signal-to-noise (S/N) problems for all but the 
brightest AGB stars. The sensitivity of Herschel is a clear im- 
provement over 750 but the short wavelength limit of PACS 
(~60 fim) is somewhat of a limitation. Nevertheless dust-species 
like Forsterite (Mg2Si04) at 69 pm, Calcite CaCOa at 92.6 pm. 
Crystalline water-ice at 61 pm, andHibonite CaAliaOig at 78 pm 
are expected to be detected. Other measured features lack an 
identification e.g. the 62-63 //m feature with candidate sub- 
stances like Dolomite, Ankerite, or Diopside (see Waters 2004 
and Henning 2010 for overviews). At longer wavelengths, PAH 
'drum-head' or 'flopping modes' have been predicted to occur 
(Joblin et al. 2002), that can be looked for with the SPIRE FTS 
(Fourier Transform Spectrometer) that will observe in an previ- 
ously unexplored wavelength regime. 

Apart from soUd state features the PACS and SPIRE range 
contain a wealth of molecular lines. Depending on chemistry 
and excitation requirements, the different molecules sample the 
conditions in different parts of a circumstellar envelope (CSE). 
While for example CO observations in the J= 7-6 line (370 pm) 
can be obtained under good weather conditions from the ground, 
this line traces gas of about 100 K. With SPIRE and PACS one 
can detect CO J= 45-44 at 58.5 pm at the short wavelength 
edge of PACS (as was demonstrated in Decin et al. 2010a) 
which probe regions very close to the star Although only the 
Heterodyne Instrument for the Far Infrared (HIFI, de Graauw 
et al. 2010) onboard Herschel will deliver resolved spectral line 
observations, PACS and SPIRE with their high throughput will 
allow full spectral inventories to be made. The analysis of PACS, 
SPIRE (and HIFI and ground-based) molecular line data with 
sophisticated radiative transfer codes (e.g. Morris et al. 1985, 
Groenewegen 1994, Decin et al. 2006, 2007) will allow quanti- 
tative statements about molecular abundances, the velocity struc- 
ture in the acceleration zone close to the star, and (variations in) 
the mass-loss rate. 

With these science themes in mind, the preparation for a 
Guaranteed Time (GT) Key Program (KP) started in 2003, 
culminating in the submission and acceptance of the MESS 
(Mass-loss of Evolved StarS) GTKP in June 2007. It involves 
PACS GT holders from Belgium, Austria and Germany, the 
SPIRE Specialist Astronomy Group 6, and contributions from 
the Herschel Science Centre, and Mission Scientists. The allo- 
cated time is about 300 hours, of which 170h are devoted to 
imaging and the remaining to spectroscopy. 

Section 2 describes the selection of the targets and Section 3 
describes the observing strategy. Section 4 discusses some as- 
pects of the current data reduction strategy. Section 5 presents 
the key science topics that will be pursued and this is il- 
lustrated by highlights of the results obtained in the Science 
Demonstration Phase (SOP), and presenting ongoing efibrts. 
Aperture photometry for 70 AGB and post-AGB stars is pre- 
sented and compared to AKARI data. Section 6 concludes this 
paper. In two appendices details on the PACS mapping and data 
reduction strategy are presented. 

2. Target selection 

2.1. AGB stars and Red SuperGiants 

The main aim of the imaging program is to resolve the CSEs 
around a representative number of AGB stars, and thereby study 
the global evolution of the mass-loss process and details on the 
structure of the CSE. With a typical AGB lifetime of lO*" year 
and a typical expansion velocity of 10 km s ' (see Habing & 
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Olofsson 2003) the effects of the mass-loss process could, in 
principle, be traced over 3 ■ lO'"* km, or 10 pc, or about 30' at 
1 kpc distance. In practice the outer size of the AGB shell will 
be smaller, first of all due to interaction of the expanding slow 
wind with the ISM, and by observational limits in terms of sen- 
sitivity and confusion noise. 

The starting point of the target selection was the ISO archive 
from which all objects classified as "stellar objects" with SWS 
and LWS observations, as well as all sources from programs 
which had "AGB stars" in the proposal keyword, were compiled. 
In addition, stars showing extended emission in the IRAS 60 or 
100 jum bands (Young et al. 1993) were considered as well. From 
that, a master list of about 300 objects was selected of stars seem- 
ingly AGB stars or related to the AGB, based on the spectral 
type, and/or simbad classification. 

The final sample was chosen to represent the various types 
of objects, in terms of spectral type (covering the M-subclasses, 
S-stars, carbon stars), variability type (L, SR, Mira), and mass- 
loss (from low to extreme) within an overall allocated budget 
of GT for this part of the program. In the selection the IRAS 
CIRR3 flag was considered to avoid regions of high background. 
Within each subclass, typically the brightest mid-IR objects were 
chosen. 

A sample of 30 O-rich AGB stars and RSGs, 9 S-stars, and 
37 C-stars will be imaged with PACS, as well as the two post- 
RSGs (IRC +10 420 and AFGL 2343). A subset of respectively, 
11, 2 and 13 AGB/RSG stars wifl be imaged with SPIRE, as 
well as R CrB, the prototype of its class (see Table [TJ. That the 
PACS and SPIRE target lists are not identical is on the one hand a 
question of sensitivity-the fluxes are expected to be higher in the 
PACS wavelength domain-and on the other hand is a question of 
the available hours of guaranteed time available to the different 
partners. 

The targets for PACS and SPIRE spectroscopy are (with one 
exception) a subset of the imaging targets. They have been se- 
lected to be bright with IRAS fluxes S'eoiSO and Su^^AQ Jy, a 
S/N of > 20 on the continuum is expected over most of the wave- 
length range (~55 to ~180 jum). A sample of 14 O-rich AGB 
stars and RSG, 3 S-stars, and 6 C-stars will be observed spectro- 
scopically with PACS, as well as the two post-RSGs. A subset of 
respectively, 5 O-rich and 4 C-rich AGB stars and 1 post-RSG 
star will be observed spectroscopically with SPIRE. 

2.2. post-AGB stars and PNe 

The aim of the imaging part is very similar to that for the AGB 
stars: mapping a few infrared bright objects in the post-AGB 
(P-AGB) and PNe phases of evolution, and to trace the flux and 
geometry of the earliest mass-loss in the phases. The objects pro- 
posed here are very well studied in broad wavelength regimes, 
but a clear understanding of the structures of their dust shells 
(and how these came about) is still lacking. 

A sample of twenty well-known C- and O-rich P-AGB and 
PNe will be imaged with PACS, and a subset of 9 with SPIRE 
(see Table[T]i- 

With respect to spectroscopic observations, the aims are 
again similar to those for AGB stars, and are related mainly to 
molecular chemistry and dust features. However, another theme 
will be the exploitation of the higher spatial resolution spec- 
tral mapping possibilities of PACS, and to spatially resolve fine- 
structure (FS) lines of hotter P-AGB stars and some PNe. When 
the central object becomes hotter than Teff of around 10 000 K, 
FS lines become apparent (e.g. Liu et al. 2001, Fong et al. 2001), 
and even at lower effective temperatures FS lines may arise from 



shocks (e.g. HoUenbach & McKee 1989). Lines in the PACS do- 
main include: [N II] 122 fim and 205 yum, [O III] 88 fim, [O I] 
146 fim, [C II] 158 fim. The first three lines come from ionised 
regions while the fluxes of the last two are emitted by photodis- 
sociation regions (PDRs). For the hotter PNe (>20-25 kK), the 
FS lines from the ionised regions are good tracers of the local 
conditions and abundances (e.g. Liu et al. 2001). For instance, 
the [N II] 122/205 fim line-ratio gives a good tracer for the elec- 
tron density, which is not very dependent on the electron tem- 
perature. FS lines emitted in PDRs provide a way to determine 
directly the PDR temperatures, densities, and gas masses. The 
spatial resolution of PACS may enable us to locate the ionised 
regions and PDRs in the nebulae. 

In total, twenty-three C- and O-rich P-AGB stars and PNe 
will be observed spectroscopically with PACS, and a subset of 
1 1 with SPIRE, of which 9 are in common with PACS. 

2.3. Massive stars 

Although massive O-type stars have intense stellar winds, their 
mass-loss rates are not sufficient for them to become WR stars 
without assuming short episodes of much stronger mass-loss 
(Fullerton et al. 2006). These short-lived evolutionary stages of 
massive star evolution (~ lO^yrs for the LBV phase) produce 
extended regions of stellar ejecta. In fact, most LBVs and many 
WR stars are nowadays known to be surrounded by ring nebu- 
lae consisting either of ejecta or ISM material swept up by the 
fast stellar winds of the stars (e.g. Hutsemekers 1994, Nota et al. 
1995, Marston 1997, Chu 2003). These stellar winds are key to 
understand the stellar evolution, as well as galactic population 
synthesis (e.g. Brinchmann 2010) 

Most of these nebulae contain dust (Hutsemekers 1997) and 
many of them were detected with IRAS at colour temperatures as 
low as 40 K. However, most of these objects were not resolved 
with IRAS, whereas the angular resolution of the instruments on 
board Herschel (6-36") will provide accurate maps of their far 
IR emission. This should help to infer detailed information on 
the mass, size and structure of the dust shells, on possible grain 
size / temperature gradients, therefore significantly improving 
our knowledge of the mass-loss history of the crucial LBV-WN 
transition phase. 

PACS photometry will be obtained for several targets among 
flie most representative LBV (AG Car, HR Cai-, HD 168625) and 
WR nebulae (M 1-67, NGC 6888). These targets represent vari- 
ous types of LBV and WR nebulae at different stages of interac- 
tion with the ISM. PACS spectroscopy will be obtained for the 
brightest nebulae to measure the fine-structure lines and deter- 
mine the physical conditions and abundances in the ionised gas 
and/or in the photodissociation region. 

In total, eight stars will be mapped with PACS, 2 will be ob- 
served with PACS spectroscopy and 2 with SPIRE spectroscopy 
(see Table[T] and details in Vamvatira-Nakou et al. 201 1). 

2.4. Supernova Remnants 

SCUBA observations at 450 and 850 ^m suggested that 0.1- 
1 Mq of cold SN dust (with a high degree of polarisation) 
was present in the Cas A (Dunne et al 2003, 2009) and Kepler 
(Morgan et al. 2003; Gomez et al. 2009) supernova remnants. No 
cold dust component was seen above the strong non-thermal syn- 
chrotron emission in SCUBA observations of the Crab remnant 
(Green, Tuffs & Popescu 2004). Arguments against such large 
quantities of cold dust being present in Cas A were advanced by 
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Krause et al. (2004) and Wilson & Batrla (2005), where they 
suggested that intervening molecular clouds along the line of 
sight to the remnant are responsible for the sub-mm emission 
seen with SCUBA, and by Sibthorpe et al. (2010). An alterna- 
tive explanation for the excess sub-mm emission at 850 yum was 
proposed by Dwek (2004), in which elongated iron needles with 
a high sub-mm emissivity could also account for the excess flux 
without the need for large dust masses. We propose to address 
this controversial and as yet, unresolved issue by searching for 
newly formed dust in five young (< 1000 yr old) Galactic SNRs 
(this should ensure they are in the late expansion or early Sedov 
phase and are therefore dominated by the SN ejecta dynamics 
and not by swept up material). We will obtain PACS and SPIRE 
far-IR and sub-millimetre photometric and spectroscopic map- 
ping of the remnants: Cas A, Kepler, Tycho, Crab and 3C 58 {- 
SN 1181); the spectroscopic insight will be important in helping 
to disentangle the issue of interstellar versus supernova and/or 
stellar wind material. 

3. Observing strategy 

3.1. SPIRE imaging 

Each of the 27 evolved star targets will be mapped in "Large 
Map" mode, with a scan leg length of 30', a cross-scan length 
of 30' and a repetition factor of 3. Such an observation takes 
4570 sec of which 2800 sec are spent on-source. The 5 SNe rem- 
nants will be mapped in the same mode and repetition factor over 
32'x32' to ensure sufficient sky coverage. This should yield 5cr 
detections for a point source with fluxes of 26, 22 and 3 1 mJy 
at 250, 350 and 500 fim, or 5cr per beam for an extended source 
with surface brightnesses of 6.0, 4.5 and 2.0 MJy sr"' at these 
wavelengths, and will reach the confusion noise at these wave- 
lengths. These numbers, and the sensitivities quoted in the next 
three sub-sections, are based on HSPOT Herschel Observation 
Planning Tool v5.1.1, which includes calibration files based on 
the in-flight performance of the instruments. 

One object (the Helix Nebula) is so large that it will be ob- 
served in SPIRE-PACS parallel mode over an area of about 75' 
squared. 

As some imaging data obtained in SDP was made public, a 
few hours of observing time could be put back into the program 
as compensation. Based on the PACS observations available in 
May 2010, it was decided to obtain a 8'x8' SPIRE map of Ml- 
67 (a massive object), a 15'xl5' map of NGC 6853 and lO'xlO' 
maps of NGC 650 and 8 AGB targets. 

The flux calibration uncertainty is 15% for all SPIRE pho- 
tometer bands at the time these data were reduced (Swinyard et 
al. 2010). The three SPIRE beams are approximately Gaussian 
with FWHM sizes of 18.1, 25.2, and 36.6", and an uncertainty 
of 5%. 

3.2. SPIRE spectroscopy 

A single pointing FTS observation will be obtained for each of 
the 23 evolved star targets, with sparse image sampling. The 
high-Hlow spectral resolution setting (the highest unapodised 
spectral resolution is 0.04 cm"') will be used, with a repetition 
factor of 17 for each mode. 

For a source with fluxes of 24, 5 and 1 .8 Jy at 200, 400 and 
600 fj.m, this is predicted to yield continuum S/N ratios per reso- 
lution element of 57, 18 and 4.8, respectively, in the high spectral 
resolution mode. Such an observation takes 2710 sec in total, of 
which 2260 sec are spent on-source. 



For the Cas A supernova remnant, FTS spectroscopy in high 
resolution mode will be obtained at three positions on the rem- 
nant, each sparsely sampling the 2.6' field of view of the FTS, 
with a repetition factor of 24. Similar spectra will be obtained 
for one pointing position only for each of the Crab and Tycho 
supernova remnants. 

3.3. PACS imaging 

All PACS imaging is done using the "scan map" Astronomical 
Observation Request (AOR) with the "medium" scan speed of 
20"/s (with the exception of the Helix nebulae mentioned ear- 
lier that will be observed using PACS-SPIRE parallel mode at 
60"/s speed). Our observations are always the concatenation of 
2 AORs (a scan and an orthogonal cross-scan). 

The P-AGB and PNe objects will be imaged at 70 and 160 
jum. The size of the maps range from about 4 to 9 arcmin on a 
side with a repetition factor of 6-8. The five SNe remnants will 
be observed with 22' scan-leg lengths at 70, 100 and 160 fim 
with an repetition factor of 2. 

For the AGB stars, the largest sub-sample within the pro- 
gram, dust radiative transfer calculations have been performed 
for all sources with the code DUSTY (Ivezic et al. 1999), to 
predict the total flux and the surface brightness as a function of 
radial distance to the star at the PACS and SPIRE wavelengths 
under the assumption of a constant mass-loss rate, taking litera- 
ture values for the stellar parameters, like effective temperature 
and distance. The current mass-loss rate and the grain properties 
were determined by fitting the SED, using archival ISO SWS 
+ LWS spectra (or IRAS LRS spectra if no ISO spectra were 
available) and broad-band photometry from the optical and infra 
red and including sub-mm fluxes when available, as constraints. 
This work is described in detail in the Ph.D. thesis by D. Ladjal 
(2010). The model predictions were then compared to the back- 
ground and confusion noise estimates coming from HSPOT as 
well as by inspecting the IRAS Sky Survey Atla£] (Wheelock et 
al. 1994) 100 fim maps to arrive at an optimal map size. In the 
end, the source size was determined by the radius where a Icr 
noise at 70 /im of 3 mJy/beam was predicted. Before submission 
of the final AORs, we also took advantage of the fact that images 
obtained with Spitzer and A£4/?/ revealed interaction of the CSE 
with the ISM for some targets (Ueta et al. 2008, 2010, Izumiura 
et al. 2009) and these set the minimum size of the maps. 

The AGB objects will be imaged at 70 and 160 fiva and the 
scan-length varies between 6 and 34'. Repetition factors between 

2 and 8 are used. For reference, a single scan over a typical area 
with a repetition factor of 1 results in a 1-cr point-source sensi- 
tivity of 12, 14, 26 mJy or 1 1.9, 1 1.5, 5.3 MJy sr"' in the 70, 100 
and 160 /im bands, and this is expected to decrease with 1/ ^|n for 
the number of repetitions considered here. As a scan and cross- 
scan are always concatenated there is another factor V2 gain in 
sensitivity. The duration of a concatenated scan plus cross-scan 
for a typical observation with a scan-leg of 15', a cross-scan step 
of 155" (see Appendix A for details) and a repetition factor of 

3 takes 2420 sec of which 1350 sec are spent on-source. Every 
point in the sky located in the central part of the map that is cov- 
ered approximately in a homogeneous way (see Appendix A for 
details) is observed theoretically for about 31 sec. The typical 
confusion noise in the three PACS bands are 0.2-0.4, 0.5-1.2, 
0.8-1.8 MJy sr"', which implies that the deepest maps we will 
carry out may just reach the confusion noise in the reddest band, 
but never at shorter wavelengths, contrary to SPIRE where the 

' http://irsa.ipac.caltech.edu/applications/IRAS/ISSA/ 
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confusion limit is reached in all bands after a few (> 3) repeti- 
tions (Nguyen et al. 2010). 

The flux calibration uncertainty currently is 10% for the 70 
and 100 fim bands and better than 15% at 160 fim (Poglitsch 
et al. 2010). The PACS beams are approximately Gaussian with 
FWHM sizes of 5.6, 6.8, and 11.4", at 70, 100 and 160 fim, re- 
spectively (at a scanspeed of 20"/s). At a low flux level a tri-lobe 
structure in the PSF is seen (see Fig lB.ll ). that is explained by the 
secondary mirror support structure, see Pilbratt et al. (2010) and 
Poglitsch et al. (2010). 

3.4. PACS spectroscopy 

All PACS spectroscopy is performed by concatenating one (or 
more) (B2A-Hshort Rl) and (B2B-Hlong Rl) AORs to cover the 
full PACS wavelength range from about 54 to 210 fim with 
Nyquist wavelength sampling. The duration of these 2 concate- 
nated AORs is 3710 sec of which 3130 sec are spent on-source. 
The observing mode is pointed with chopping/nodding with a 
medium chopper throw of 3'. 

The sensitivity of the PACS spectrometer is a complicated 
function of wavelength (see Section 4 in the PACS User Manual, 
which is available through the Herschel Science Centre) but is 
best in the 110-140 //m range at about 0.8 Jy (1-cr) continuum 
sensitivity for a repetition factor of one. The spectral resolution 
also strongly depends on wavelength (see Section 4 in the PACS 
User Manual) and is lowest at about R- 1 100 near 110 jim and 
highest at about R- 5000 near 70 yum. 

Four sources were submitted as SDP targets for PACS 
spectroscopy, but three (NGC 7027, VY CMa and CW Leo) 
happened to be observed already as part of the Performance 
Verification (PV) of the AORs. These data were oflicially "trans- 
ferred" to the MESS program. Two of these sources (VY CMa 
and CW Leo) were observed in mapping mode. All details over 
these maps can be found in Royer et al. (2010) and Decin et al. 
(2010b) respectively. 

Nine positions on the Cas A supernova remnant will be ob- 
served with the PACS IFU, obtaining full spectral coverage using 
the (B2A-Hshort Rl) and (B2B-Hlong Rl) AORs and the largest 
chopper throw of 6'. Single positions on each of the Crab, Kepler 
and Tycho remnants will be observed in the same way. 



4. Data reduction 

4.1. PACS imaging 

Poglitsch et al. (2010) describe the standard data reduction 
scheme for scan maps from so-called "Level 0" (raw data) to 
"Level 2" products. Level 1 and 2 products that are part of the 
data sets retrieved from the Herschel Science Archive (HSA) 
have been produced by execution of the standard pipeline steps, 
as described by Poglitsch et al. (2010) 

We do not use the products generated in this way, but use 
an adapted and extended version of the pipeline script suited to 
our needs starting from the raw data. In particular the deglitch- 
ing step(s), and the "high pass filtering" need special care. 
Therefore, the making of the map has become a two-step pro- 
cess. In addition, as the standard pipeline script operates on a 
single AOR, while our observations are always the concatena- 
tion of 2 AORs (a scan and an orthogonal cross-scan) an addi- 
tional step is also needed. Details of the reduction steps can be 
found in Appendix B. 



4.2. PACS spectroscopy 

Poglitsch et al. (2010) describe the standard data reduction 
scheme for chopped spectroscopic measurements and the spec- 
troscopic results described in Sect. |5]have been obtained using 
the nominal pipeline except that we combined the spectra ob- 
tained at the two nod positions after rebinning only (Royer et al. 
2010). 

4.3. SPIRE imaging 

The standard SPIRE photometer data processing pipeline, de- 
scribed by Griflin et al. (2008, 2010) and Dowell et al. (2010), 
is sufficient to reduce our SPIRE photometric imaging data. The 
calibration steps for these data are described by Swinyard et al. 
(2010). The high stability and low I// noise knee frequency of 
the SPIRE detectors mean that no filtering of the signal time- 
lines is required. The larger beam size and characteristic re- 
sponse of the SPIRE bolometers to cosmic rays means that the 
standard deglitching routine is applicable in this case. 

The only custom step required in the reduction of these data 
is the removal of the bolometer base-lines. The nominal process- 
ing removes only the median from each time-line. In the pres- 
ence of large scale structure and/or galactic cirrus, this method 
can result in striping due to sub-optimal baseline estimation. To 
overcome this we implement the iterative baseline removal de- 
scribed by Bendo et al. (2010). 

This method starts by creating an initial map from the me- 
dian subtracted signal time-Unes. The signal in a given bolome- 
ter time-line is then compared to the value of the derived map 
pixel for each pointing within a scan-leg. The median difference 
is then calculated, and subtracted from the time-line. This pro- 
cess is repeated 40 times to fully remove the striping. Finally the 
median is removed from the last iteration to give the final map 
product. 

4.4. SPIRE spectroscopy 

Griffin et al. (2010) and Fulton et al. (2010) describe the SPIRE 
spectrometer data processing pipeline, while Swinyard et al. 
(2010) describe how the spectra are calibrated. We briefly re- 
view here the data reduction of the MESS targets. 

Data reduction was carried out in HIPE (Herschel 
Interactrive Processing Environment, Ott et al. 2010). An impor- 
tant step in SPIRE spectral data reduction is accurate subtraction 
of the background emission from the telescope and instrument. 
During Herschel's SDP, a reference interferogram was obtained 
on the same Operational Day (OD) as target observations. The 
sources observed at these times were also relatively bright. The 
subtraction of a reference interferogram obtained close in time 
to the target spectrum gave good results. 

Since the beginning of routine observations, fainter targets 
have been observed, and deep reference interferograms have 
been obtained less frequently. If the temperature of the telescope 
and instrument differ significantly between the dark observation 
and the source observation, unwanted spectral artefacts can ap- 
pear, the most obvious being a ski-slope-like shape to the con- 
tinuum in the SLW band. 

Routine phase calibrations now use a "superdark" created by 
combining a large number of dark observations to obtain a very 
high signal-to-noise ratio. A short dark observation obtained on 
the same OD as the source spectrum is then used to scale the su- 
perdark to correct for the variations in instrument and telescope 
temperatures. 
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Data in the SDP were flux-calibrated using observations of 
the asteroid Vesta, due to the inaccessibility at that time of 
Uranus and Neptune, the preferred calibration sources. Routine 
phase flux calibration uses observations of Uranus and the plan- 
etary atmosphere models of Moreno (1998, 2010). The Vesta- 
based overall calibration accuracy was estimated at 15-30% 
(Griflin et al. 2010); the Uranus-based calibration is a significant 
improvement, particularly in the SLW band. 

The FTS produces spectra in which the line profiles are sine 
functions; to remove the negative side-lobes of the sine func- 
tions, we apodize our spectra using the extended Norton-Beer 
1.5 apodizing function described in Naylor & Tahic (2007). This 
gives line profiles which are close to Gaussian, with a full width 
at half-maximum (FWHM) of 0.070 cm ' (2. 1 GHz). 

5. Key science and first results 

The results obtained during SDP and the early phases of Routine 
Phase (RP) match and exceed our expectations. Eight papers 
were published in the A&A special issue and one in Nature that 
illustrate well the science that will be pursued. 

5.1. Supernova dust 

Barlow et al. (2010) analyse PACS and SPIRE images of the 
CC SN remnant Cas A to resolve for the first time a cool dust 
component, with an estimated mass of 0.075 Mq, significantly 
lower than previously estimated for this remnant from ground- 
based sub-millimeter observations (see Section 2.4) but higher 
than the dust masses derived for nearby extragalactic supernovae 
from Sp/fzer observations at mid-IR wavelengths (see Section 1). 

5.2. The circumstellar envelopes 

Van Hoof et al. (2010) analyse PACS and SPIRE images of the 
planetary nebula NGC 6720 (the Ring nebula). There is a strik- 
ing resemblance between the dust distribution and H2 emission 
(from ground based data), which appears to be observational evi- 
dence that H2 has been formed on grain surfaces. They conclude 
that the most plausible scenario is that the H2 resides in high 
density knots which were formed after recombination of the gas 
started when the central star entered the cooling track. 

Another main result is on the detached shells around AGB 
stars that were first revealed by the IRAS satellite as objects that 
showed an excess at 60 /^m. Later observations by single-dish 
(Olofsson et al. 1996) and interferometric millimetre telescopes 
(Lindqvist et al. 1999, Olofsson et al. 2000) revealed spherical 
and thin shefls emitting in CO (TT Cyg, U Ant, S Set, R Scl, U 
Cam), that were interpreted as short phases of high mass loss, 
probably related to thermal pulses. Kerschbaum et al. (2010) 
present for the first time resolved (PACS) images of the dust 
shells around AQ And, U Ant, and TT Cyg, which allows the 
derivation of the dust temperature at the inner radius. Although 
U Ant and TT Cyg are classical examples of stars surrounded 
by a detached shell, AQ And was not previously known to have 
such a shell. Additional examples are presented in Kerschbaum 
et al. (201 1) and Mecina et al. (201 1). 

5.3. Wind-ISM interaction 

Sahai & Chronopoulos (2010) present ultraviolet GALEX im- 
ages of IRC -1-10 216 (CW Leo) that revealed for the first time a 
bow shock. Ladjal et al. (2010) demonstrate that the bow shock 
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Fig. 1. The map of X Her at 70 fim displayed on a logarithmic 
intensity scale. North is to the top, east to the left. To the south 
east of X Her are a pair of galaxies (UGCI0156a and b). 



is also visible in PACS and SPIRE images. The dust associated 
with the bow shock has a temperature estimated at 25 K. Using 
the shape of the shock and a published proper motion, a space 
motion of CW Leo relative to the ISM of 107/ V«ism km s"' is 
derived. A comparison to the models by Wareing et al. (2007) 
regarding the shape of the bow shock suggests that the space 
motion is likely to be <75 km s"', implying a local ISM density 
>2 cm""*. In fact, examples of bow shocks are quite common in 
the MESS sample. Figure [T] shows the case of X Her (Jorissen 
et al., in prep.) which displays extended emission over about 2'. 
Young et al. (1993) quote an extension in the IRAS 60 //m band 
of more than 12' but we find no evidence for that. Interestingly, 
panel B in their Figure 1 1 shows an intensity profile for X Her 
which appears to be much narrower than the size quoted in their 
Table 3. Additionally, IRAS did not have the spatial resolution to 
resolve the pair of interacting galaxies we find, possibly leading 
to an overestimate of the size of the extended emission by Young 
et al. In the case of X Her the flattened shape of the bow shock 
suggests a larger space motion and/or a larger ISM density com- 
pared to CW Leo, or to a highly inclined orientation of the bow 
shock surface with respect to observers. Additional examples are 
presented in Jorissen et al. (201 1) and Mayer et al. (201 1). 

5.4. Molecular lines and chemistry 

Regarding spectroscopy, the results published so far illustrate the 
enormous potential to detect new molecular lines thanks to the 
increased spectral resolution of PACS w.rt. ISOs LWS and the 
previously unexplored wavelength region covered by the SPIRE 
FTS. 

Three papers concentrate on CW Leo. Cernicharo et al. 
(2010) discuss the detection of HCI in this object, and derive 
an abundance relative to H2 of 5x10 and conclude that HCI is 
produced in the innermost layers of the circumstellar envelope 
and extends until the molecule is photodissociated by interstel- 
lar UV radiation at about 4.5 stellar radii. Decin et al. (2010a) 



M.A.T. Groenewegen et al.: MESS (Mass-loss of Evolved StarS), aHerschel Key Program 



7 



discover tens of lines of SiS and SiO including very high transi- 
tions that trace the dust formation zone. A comparison to chemi- 
cal thermo dynamical equilibrium models puts constrains on the 
fraction of SiS and SiO that are involved in the dust formation 
process. The discovery of water a few years ago by Melnick et al. 
(2001) spurred a lot of interest. However the exact origin could 
so far not be identified. Decin et al. (2010b) analyse SPIRE and 
PACS spectroscopy and present tens of water lines of both low 
and high excitation. A possible explanation is the penetration of 
interstellar ultraviolet photons deep into a clumpy CSE initiating 
an active photo chemistry in the inner envelope. 

That water will be a main theme is also illustrated in Royer 
et al. (2010) who present PACS and SPIRE spectroscopy of the 
RSG VY CMa. Nine hundred lines are identified, of which half 
are of water. Finally, Wesson et al. (2010b) present SPIRE FTS 
spectra of three archetypal carbon-rich P-AGB objects, AFGL 
618, AFGL 2688 and NGC 7027, with many emission lines de- 
tected. Results include the first detection of water in AFGL 2688 
and the detection of the fundamental J- 1-0 line of CH^ in the 
spectrum of NGC 7027 (see Fig. |2]i. Figure |2] also shows the 
PACS spectrum and for comparison the ISO LWS spectrum in 
Liu et al. (1996) scaled by a factor of 0.25 in flux. The inset more 
clearly illustrates the improved sensitivity and spectral resolu- 
tion over a smaller wavelength region. The width of the LWS res- 
olution element was 0.29 or 0.6 /im (depending on wavelength), 
while that of PACS varies between 0.03 and 0.12 ^m. The LWS 
effective aperture was 80" diameter, while the PACS spectrum 
is based on the central 3x3 spatial pixels, corresponding to 28 
X 28". The continuum shape of the PACS spectrum may still 
affected by the fact the ground-based RSRF (relative spectral re- 
sponse functions) was used. The noise appears to vary across 
the spectrum because the red and blue detectors have different 
characteristics. A full analysis of the PACS spectrum, including 
a detailed comparison to the LWS spectrum will be presented 
elsewhere (Exter et al. 201 1). 

5.5. Dust Spectroscopy 

De Vries et al. (2011) will present the first result on dust spec- 
troscopy, the detection of the 69 fim Forsterite feature in the P- 
AGB star HD 161796. Identifying dust features has proven to 
be perhaps not as easy as originally believed. Reasons are that 
a more accurate determination of the PACS and SPIRE RSRFs 
and improved data reduction techniques are still ongoing, and 
also the wealth of molecular and fine-structure lines that need to 
be removed before relatively broad dust features can be identi- 
fied. 



5.6. AGB and post-AGB stellar fluxes and a comparison to 
literature data 

The imaging observations of the AGB and P-AGB stars in the 
MESS program have been largely completed for SPIRE and 65 
percent completed for PACS, and in this section we present the 
fluxes of the central objects. The fluxes are calculated by aper- 
ture photometry. The size of the aperture is calculated in an auto- 
mated way in a similar way as the determination of the source ra- 
dius described in Appendix B: The rms in the final map is deter- 
mined, and then the 15, 12, 8, 12, 10, 8 cr contours (for PACS 
70, 100, 160, SPIRE 250, 350, 500 /urn, respectively) are approx- 
imated by a circle. This radius of this circle can be multiplied by 
a user-supplied value. The default value is 1 .4 so that the aper- 
ture includes all emission above approximately 3cr. For objects 
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Fig. 2. Top panel: SPIRE FTS spectrum of NGC 7027 (Wesson 
et al. 2010b), with '-CO, '^CO and CH+ Unes indicated. Bottom 
Panel: PACS spectrum of NGC 7027 (Exter et al. 201 1 ; the blue 
upper curve) with the strongest lines indicated, and the ISO LWS 
spectrum (black, lower curve) of Liu et al. (1996) scaled by a 
factor of 0.25 in flux. The inset shows a detail of the spectrum 
(with no scaling in flux) illustrating the improved resolution of 
the PACS spectrum w.rt. the LWS spectrum. 



with detached shells or extended emission due to wind-ISM in- 
teraction the aperture was changed manually to include as best as 
possible the central object only (see the remarks in Table|2]i. The 
location of the sky annulus for subtraction of the background 
was verified manually to be as much as possible free of extended 
emission. 

We tuned our automatic procedure of choosing the apertures 
on PACS data of the asteroid Vesta taken on OD 160 during PV 
and the empirical beam profiles provided by the SPIRE team. We 
find apertures of 28 (PACS 70 fim) and 36 (PACS 100 fim) arc- 
second for which the PACS Scan Map Release Note states that 
they enclose about 97% of energy. For the three SPIRE bands 
we find radii of 64, 80 and 118". The ratio of this size to the PSF 
size is similar to that for PACS, suggesting that these radii also 
include about 97% of energy. 

Table|2]lists the fluxes of the central source for the AGB and 
P-AGB stars observed by October 17th. Objects where the fluxes 
of the central object may still include emission from extended 
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emission are flagged. Listed are the flux and the aperture used 
between parenthesis. The formal error on the flux is negligible 
and almost always below 0.01 Jy (PACS), and 0.02 Jy (SPIRE). 
The absolute flux calibration accuracy of the PACS and SPIRE 
photometry is estimated to be 10-15% (Poglitsch et al. 2010, 
Griflin et al. 2010). 

To make a first consistency check of the derived photometry 
we compare the PACS 70 and 160 //m data with each other. The 
mean flux ratio between 70 and 160 jjm of the observed sources 
is 6.1, with standard deviation of 1.5, which is a little higher, but 
within the standard uncertainty, than the value of 5.2 expected 
for the ratio in fluxes of a modified black body (v^By(T), with 
j8 = 2) on the Rayleigh- Jeans tail at the observed wavelengths 
of 70 and 160 fim. However, the scatter is relatively large and 
for 13 out of the 68 sources the flux ratio deviates more than 
30% (Icr) from the mean. The targets which have a significantly 
lower flux ratio are all carbon-rich AGB or P-AGB, while the 
stars with higher ratios are all oxygen-rich AGB (and 1 oxygen 
rich P-AGB). If we split the oxygen-rich and carbon rich stars 
in two groups, the mean 70/160 flux ratios are 7.3 +1.4 and 
5.5 ± 1.2, respectively. Although both are consistent with the 
expected value, these ratios could indicate a larger value of /3 
for the oxygen-rich stars than for the carbon-rich stars. These 
effects are illustrated in Fig. [3] On the other hand a systematic 
difference in the dust temperature could alter these ratios as well, 
in particular if the spectral energy distribution deviates from the 
Rayleigh- Jeans approximation at these wavelengths. 

A second consistency and data quality check of the mea- 
sured PACS central source ffux densities is performed by cross- 
matching targets and comparing ffuxes in this work with the 
AKARI/FIS bright source catalog (Yamamura et al. 2010). We 
find 59 matches with the FIS bright source catalog, of which 49 
have good quality 65 //m data (the low quality data correspond 
to FIS sources brighter than 500 Jy or fainter than ~2 Jy). The 
correlation coefficient for these two data sets (omitting the bad 
quality FIS flux densities) is very good with R- 0.98 (see Fig.|4|. 
The regression improves (in particular at lower flux levels) from 
choosing a power law instead of a linear function. The slope of 
the linear function is close to expected flux ratio of 1 . 1 (again for 
a modified black body with /3 - 2). NML Tau and o Cet deviate 
significantly from this relation but these sources are relatively 
bright (>200 Jy) and may therefore suffer from inaccurate FIS 
fluxes. Other sources deviating from the linear relationship (with 
FIS fluxes higher than PACS fluxes) are NML Tau, WX Psc, and 
U Ant (see below). A similar comparison of PACS 70 fim flux 
densities with IRAS 60 jum values is shown in Fig. |5] We note 
that in general flux densities appear consistent with previous re- 
sults thus warranting a direct comparison for SED modelling. 
However in a few cases extended or detached emission near the 
AGB stars contaminated previous flux measurements of the cen- 
tral source (see also below for examples), thus introducing errors 
in the results obtained with SED models. Additional errors could 
arise from flux variability on time-scales of a few years up to 
two decades. Also, errors in the automated point source extrac- 
tion routines for both IRAS and AKARI could have introduced 
unknown errors for particular objects. 

A comparison between the PACS 160 jim and FIS 140 or 
160 yum flux densities proved less conclusive due to the relative 
large uncertainties on the FIS fluxes (these are not reliable for 
reported fluxes below ~3 Jy, nor above 100 Jy). Provisionally, 
the PACS 160 yum fluxes are consistent with previous measure- 
ments though in several cases we noticed that the measured val- 
ues are higher than expected from the overall SED shape. Again, 
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Fig. 3. PACS 70 fim versus PACS 170 jum flux densities for the 
observed targets in Table 2. 



the presence of strong emission lines could affect the broad band 
flux measurements. 

We discuss a few particular cases from the comparisons 
made above. For the S-type star RZ Sgr the 70 fim flux is only 
1.5 times lower than at 160 /im. It is likely that the latter value 
is still contaminated by extended emission in the aperture. The 
70 ;um flux (4.2 ± 0.4 Jy) is consistent with the 65 |/m FIS flux 
(4.7 ±1.1 Jy). Also, the FIS fluxes at 140 and 160 fim are some- 
what higher, though of low quality, than the PACS 160 //m flux 
(i.e. ~4 - 5 Jy versus ~2.9 Jy). Likely, the FIS aperture is affected 
by the extended emission seen in the PACS red image. TT Cyg 
also has an average 70/160 /im flux ratio. The central source flux 
at 70 i^m reported for PACS is a third of the FIS 65 /im flux. Here 
the PACS 70 //m map shows a bright detached shell. Therefore 
the FIS flux could include both the central source and detached 
shell. Additionally, the FIS 65 //m flux has a low-quality flag. 
In this case the FIS fluxes at 140 and 160 fim (i.e. 1.5 and 0.3 
Jy) agree better with the 160 /im flux of 0.14 Jy. The detached 
shell is fainter at 160 fjm and thus the difference between PACS 
and FIS will be less. The PACS 70 and 160 fim point source 
fluxes for U Ant (6.3 and 2.5 Jy, respectively) are significantly 
lower than the corresponding FIS fluxes of 14.0, 14.0, 5.8 and 
4.8 Jy at 65, 90, 140, and 160 fim. In this case we surmise that 
the bright extended structure visible in the PACS images (which 
is not included in the point source flux) is included in the point 
source flux reported in both the FIS and IRAS catalogs. For AQ 
And, which has a bright detached shell in both PACS images, 
the central source is a factor two to three fainter than the cor- 
responding fluxes reported in the FIS catalog. However, for this 
faint source the FIS values have a low quality flag. Nevertheless, 
it seems likely that the latter include (at least part of) the shell 
emission as resolved by PACS. These examples are primarily 
given to illustrate the complexities involved when comparing re- 
sults from different surveys/instruments and in particular when 
those results are used to construct and model spectral energy dis- 
tributions. The latter will benefit also from additional informa- 
tion provided by previous missions such as ISO and Spitzer, as 
well as PACS spectroscopy observations of some MESS targets 
(Table 1). 
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Fig. 4. AKARI/FIS 65 nva flux density plotted against PACS 
70 /im flux density. The blue dashed line shows the best lin- 
ear regression, after 3cr outlier (blue open squares) rejection 
iy - 13 + l.l ■ x), while the orange solid lines is the best fit 
using a power law relation (y - 1.3 jc"'^^). The former fits were 
not constrained to go through the origin. The dotted green line 
represents the nominal flux ratio of an AO star for effective wave- 
lengths of 62.8 and 68.8 for FIS and PACS, respectively. Objects 
indicated with red filled circles were omitted from the fit due to 
low quality flag in the FIS catalog. 




PACS 70 micron flux [Jy] 

Fig. 5. PACS 70 jum versus IRAS 60 /zm flux densities for the 62 
observed targets in Table 2 with counterparts in the IRAS point 
source catalog. The solid line shows the predicted flux ratio of 
1.3 (see text) and the dashed line the powerlaw fit to all objects 
brighter than 10 Jy. 



6. Conclusions and outlook 

The scope, aims and status of the Herschel Guaranteed Time 
Key Program MESS (Mass-loss of Evolved StarS) are presented. 
The current concepts on the data processing are presented, and 
aperture photometry of all 70 AGB and post-AGB stars observed 
as per October 17th 2010 are presentecQ. 



^ The progress of the MESS project can be followed via 
WW w.univie . ac . at/ space/MES S 



Some of our SDP data is akeady public (see the remarks in 
Table [U and the data taken in routine phase have a proprietary 
period of one year, implying that data will become successively 
public from about December 2010 onwards. 

Currently, the PACS maps are produced via the standard 
PhotProject task, that use data that are filtered to remove 1/f- 
noise as input. We achieve good results for compact objects, but 
unmasked large-scale structures in the background are aff'ected 
by the filtering. In order to improve this, we are currently inves- 
tigating other filtering and mapping methods, such as MADmap 
(Cantalupo et al. 2010). That way we also seek to improve the 
spatial resolution of the final maps. A second point under in- 
vestigation is to correct the effects caused by the instrument 
PSF. With its tri-lobe pattern and other wide-stretched features 
it is currently not possible to make any definite statements about 
structure in the circumstellar emission close to the central object, 
although many sources are extended. Thus we are investigating 
different deconvolution strategies and PSF-related matters. The 
efforts the MESS consortium are currently undertaking to im- 
prove the PACS data processing are described in Ottensamer et 
al. (2011). 

On the science side, the publication of the very first Nature 
paper based on Herschel results by Decin et al. (2010b) is a 
highlight. It demonstrates the power of the PACS and SPIRE 
spectrometers. With more than 50 PACS and almost 30 SPIRE 
targets to be observed spectroscopically this will result in an ex- 
tremely rich database that, with proper modelling, will allow de- 
tailed studies on molecular abundances, the velocity structure in 
the acceleration zone close to the star, and the mass-loss rate. 

On the imaging side the fact that bow shock cases are ubiqui- 
tous is extremely interesting. Although this in fact makes it more 
difficult to derive the mass-loss rate history of the AGB star, it 
offers an unique opportunity to use these cases as probes of the 
ISM. 
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Table A.l. Scan map observing mode parameters 



Fig. A.l. Top: Coverage map of the scan and cross-scan (left and 
right) of TT Cyg obtained in SDP on OD 148. Bottom: Coverage 
map of the scan and cross-scan of AQ And obtained in Routine 
Phase on OD 224. Until OD 221 the coverage of the scan map 
obtained at an angle of 45° was not uniform. 



Appendix A: Lessons from the SDP 

A few sources were proposed to be observed during Science Demonstration 
Phase (SDP). This was intended as a demonstration of the implementation of 
the concatenation of a scan and a cross-scan, at the recommended array-to-map 
angles of 45 and 135°. It was discovered that the coverage of maps obtained at 
45° was not uniform. In the example shown in top left panel of Fig. lA. li the 
coverage in the center of the map is only about 80% of the maximum coverage. 
This was reported to the instrument team. The plate scale of the PACS camera 
was then recalibrated, which led to a 3% larger pixel size on the sky. An extra 
rotation of 2.5 degrees of the instrument coordinate system w.rt. the satellite ref- 
erence also had to be introduced. The changes have been implemented on OD 
22 1 , and since then the coverage is now much more uniform (see lower panel of 
Fig. lA. l\ . No changes are visible to the user and so in HSPOT angles of 45 and 
135 degrees still need to be entered. 

Oiiginally we opted to use "homogeneous coverage" and "square inap" as 
observing mode parameters for the scan maps. In that case the HSPOT internal 
logic takes the scan-leg length from the user input to compute the cross-scan step 
and the number of scan-legs needed. It turned out however that a small change 
in scan-leg length can then make a large difference in the computed observing 
time. Therefore, we have computed the discrete set of optimum solutions for the 
scan-leg length and cross-scan distance con'esponding to the square maps ob- 
tained from perpendicular scans performed under array-to-map angles of 45 and 
135 degrees. The standard cross-scan distance for uniform coverage is 155 arc- 
sec. For repetitions of 2 or 4, this can be divided accordingly, and intermediate 
solutions are found (see Table lA. H . which are more efficient than repeating the 
whole scan. For example, if one would like to cover an area of 60'x60' homo- 
geneously, one can opt for a cross-scan step of 155", a scan-leg length of 64.4', 
and 25 scan legs. 



Appendix B: PACS image reduction 

In this Appendix the strategy that was used to reduce the PACS images is outlined 
in some detail. It is mainly intended for persons who are interested in, or want to 
learn more about, PACS image reduction using the "high pass filter" technique. 

As mentioned in Sect. 14.11 the data reduction starts from "Level 0" (raw 
data) as retrieved from the Herschel Science Archive. Compared to the "Level 
2" product as produced by execution of the standard pipehne steps, as described 
by Poglitsch et al. (2010), special care need to be given to the deglitching step(s), 
and the "high pass filtering". Therefore, the making of the map has become a 
two-step process. In addition, as the standard pipeline script operates on a single 
AOR, while our observations are always the concatenation of 2 AORs (a scan 
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and an orthogonal cross-scan) an additional step is also needed. All steps below 
are valid and tested under HIPE v4.4.(0(Ott et al. 2010). 

The generation of the level 1 product follow the description in Poglitsch et 
al. In particular, the deglitching is done using the photMMTDeglitching task0 
and frames detected as glitches are only masked without trying to interpolate the 
signal. When the central source is bright however, then this task also incorrectly 
masks a significant number of frames of the central source. Therefore a second 
step is required as detailed below. 

In the case of a scan map AOR the level 2 product is a map. The processing 
steps from the level 1 product to this map are the highpassFilter task and the 
photProject task. 

The purpose of the high pass filter is to remove detector drift and 1/f noise. 
At the moment the task is using a median filter which subtracts a running me- 
dian from each readout. The filter box size can be set by the user For a certain 
high pass filter width (HPFW), the median is run over (2-HPFW+l) consecu- 
tive readouts, and only uses readouts that are not masked. The filter works in the 
time domain. However this corresponds to filtering on a certain spatial scale as 
the scan speed (20"/s in our case) and the sampling of the photometer signal (10 
Hz) are known. In the making of the first map a HPFW of half of the scan-leg of 
the map is used. 

The map making task simply reprojects the R.A. and Declination of every 
frame onto the sky using a certain pixel size. For the generation of this first map 
the native pixel sizes are used (3.2" in the blue, and 6.4" in the red). 

The map generated in this way (and essentially equivalent to the level 2 
product retrieved from the HSA) is then analysed in order to make the second 
and final map for an AOR. In particular the location of all significant emission 
in the map will be determined, especially the location and extent of the central 
source will be characterised. With this information the second level deglitching 
can be performed, the significant emission in the map will be masked, and the 
optimised value of the HPFW will be determined. 

The importance of properly masking the significant emission in the map and 
the choice of the HPFW (a related issue) will be illustrated in detail first. For this 
we use the observations of AQ And at 70/im. This star has a very nice detached 
shell and the PACS observations are described in Kerschbaum et al. (2010). 

In earlier versions of HIPE (vl .2) a HPFW as low as 20 was used as default. 
The final map (adding scan and cross-scan) with such a low value and no mask- 
ing is displayed in Fig. IB. II The eft'ect of "shadowing" that is seen in many level 
2 product maps where there is a bright source in the middle is clearly visible. 
The reason is obvious: without any precaution the median filter will subtract too 
much when it passes over locations of significant emission, in particular a bright 
central source. Therefore one sees dark stripes beside the central source in the 
two scan directions, as well as around the detached shell. 

The azimuthally averaged ID-intensity distribution is shown in Fig. lB.2l as 
the black solid line, illustrating in a different way the regions of negative flux in 
this map. 

The solution to this problem is to mask the central object. The task 
photReadMaskFromlmage can be used to mask all frames where the signal ex- 
ceeds a user given value. In earlier versions of HIPE however, this threshold was 
not a free parameter but fixed at a large value of 0.5 Jy/pixel. Therefore an alter- 
native procedure was developed early on. In a first step, the rms, o", in the map 
is determined in a circle of 30" diameter located | of the scan-leg length to the 
North of the map centre, 200" in this particular case. To exclude any possible 
true sources and undetected glitches in that area the rms is calculated by first 
eliminating the lowest and highest 2% of the data values. 

In a second step the contour at a level of Vcr is determined. This is on purpose 
a large value to ensure a relatively well behaved shape of the contour The longest 
contour near the centre of the map is then approximated by a circle. The centre 
of this circle is taken as the centre of gravity of the contoujj The distances from 
the centre to all elements of the contour are calculated and the radius is taken as 
the mode (hereafter sourceradius.). 

This radius can then be increased by a user-supplied value. The default value 
is 1.3, which means for most sources that a region around the central star down 
to a level of approximately 3o" is masked. 

In the case of AQ And, the region that is masked has a radius of 9.8". The 
effect is shown in Fig. lB.2l as the red dashed line, for a HPFW of 20, 100 and 150 
(bottom to top). It is obvious that masking only the central object with a small 
HPFW does not lead to a rehable result. Only for a HPFW of 150 a stable ID 



^ Data presented in this paper were analysed using HIPE, a 
joint development by the Herschel Science Ground Segment 
Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HlFl, PACS and SPIRE consortia. (See 
http://herschel.esac.esa.int/DpHipeContributors.shtml I 

using photMMTDeglitching(frames, scales=3, nsigma=9, 
incr_fact=2, mmt_mode='multiply', onlyMask= True) 

' See e.g. htt p://iocaLwasp.uwa.edu.au/"^pibourke7geometry7polyarea7^ 
for an algorithm. 



intensity distribution is reached. This result is not that surprising as the emission 
in the shell has not been masked. The green dashed-dotted line and the blue dot- 
ted line show the results if a region of radius 60, respectively, 80" is masked. In 
that case it is not necessary to use a very large HPFW. In fact, the only concern 
is that it is large enough that the inedian filter contains unmasked frame!0 In the 
cases shown the HPFW values are 33 (dashed-dotted line) and 43 (dotted line). 
The dift'erent intensity profiles show that reliable results (i.e. differences in data 
reduction strategies that lead to intensities that are within the photometric cali- 
bration uncertainly of 10%) can be achieved if areas in the map that coirespond 
to significant levels of emission are masked. 

In the data reduction scheme employed in the present paper the masked 
region is where the signal is above 3o" in the entire map (using the 
photReadMaskFromlmage task), and the procedure outlined above to derive a 
sourceradius. The scaling factor of this radius is 1.3 by default to roughly mask 
a region with a signal above 3(T around the central objecQ 

The sourceradius is also used in defining the details of the 2nd level 
deglitching task. As mentioned above, the photMMTDeglitching task incor- 
rectly inasks a significant number of frames due to the brightness and compact- 
ness of the central source. In a first step all frames that are masked by the phot- 
MMTDeglitching task inside a radius O.Sxsourceradius are unmasked. 

Then the IlndLevelDeglitchTask is used on the map, using 6(T-clipping 
over a running box of 41 elements. The value for the sigma-clipping and the 
size of the box were determined by comparing the glitch rate in an off-source 
area (in fact the same area as where the rms in the inap is calculated) to the 
glitch rate in the area around the source for a large number of objects. As the 
radiation events should be distributed randomly the photMMTDeglitching and 
IlndLevelDeglitch tasks should roughly flag the same percentage of frames. 

After the masking of regions of significant emission and the masking of 
glitches on source the highpassFilter task is run again with the optimised 
value as outlined earlier The fraines objects for scan and cross-scan are joinecQ 
and the final map is created using a pixelsize of 1" (blue filter) and 2" (red filter). 



in particular it is calculated from HPFW= max(20, int(1.01x 
inaskfracx sourceradius/ scanspeed/ 0.1) +5)), where maskfrac is a 
user-supplied number with default value 1.3. This ensures that even in 
the worst case the median filter uses 1 1 data values to calculate the me- 
dian. The minimum value of 20 units is enforced for cases where the 
central source is weak and no or no reliable sourceradius can be deter- 
mined 

' As clearly demonstrated above, the scaling factor must be larger for 
the stars with detached shells or obvious regions of extended emission 
in order to compute the appropriate value of the HPFW and maskfrac 
factors of 2.0 (U Ant, EP Aqr, o Cet, n Gru), 2.5 (TX Psc, X Her)', 3.0 
(9 Aps, T Mic, RT Vir), 3.5 (a Ori) 4.0 (R Hya), 5.0 (U Cam, S Cep, 
RZ Sgr), 6.0 (AQ And, TT Cyg), 8.0 (UX Dra) were used. 

using the join method on a Frames object: frames.join(framesl) 
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(c) PACS 201 
MESS CTKP 



Fig. B.l. Top panel: The final map of AQ And at 70 jim when 
the central star and shell are not masked during median filtering 
with a "high pass filter width" of 20. The "shadowing effect" is 
clearly visible. The central star does not look round but that is 
due to the secondary mirror support structure which introduces 
a tri-lobe shape of the PSF at a low flux level. Lower panel: The 
map when the central star and shell are masked and a high pass 
filter width of 35 is used. 
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Fig. B.2. Azimuthally averaged intensity distribution for AQ 
And at 70 /im when the central star and shell is not masked be- 
fore the median filtering, and using a "high pass filter width" 
(HPFW) of 20 (black solid line), The red dashed lines indicate 
the result if a region of 9.8" radius around the central star is 
masked using a HPFW of 20, 100 and 150 (bottom to top). The 
green dashed-dotted line and the blue dotted line show the re- 
sults if a region of radius 60", respectively, 80" is masked, with 
a HPFW of 33, respectively 43. 
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Table 1. Target List. The spectral type comes from the SIMBAD database. Sources are listed alphabetically per sub-class. 



IRAS name 


Identifier 


type 


PACS 


PACS 


SPIRE 


SPIRE Remarks 








Imag. 


Spec. 


Imag. 


Spec. 


O-rich AGB 


01556-H4511 


AFGL 278 


M7ni 


V 








17411-3154 


AFGL 5379 


M 




V 






05524+0723 


a ORl / Betelgeuse 


Mllab: 


V 


V 


V 


V 


21439-0226 


EP AQR 


M8nivar 


V 








20077-0625 


IRC -10 529 /AFGL 2514 


M: 


V 


V 






21419+5832 


juCEP 


M21ae 


V 




V 






NMLCYG 


M6ine 


V 


V 






03507+1115 


NML TAU / IK Tau 


M6me 


V 


V 


V 


V 


02168-0312 


oCET 


M7ine 


V 


V 


V 


V 


01304+6211 


OH 127.8+0.0 


M 


V 








18348-0526 


OH 26.5+0.6 


M 


V 


V 






23412-1533 


R AQR 


M7111pevar 


V 




V 




23558+5106 


RCAS 


M7111e 


V 


V 


V 




10580-1803 


RCRT 


M7111 


V 








04361-6210 


RDOR 


M8ine 


V 


V 


V 




13269-2301 


RHYA 


M7ine 


V 




V 




09448+1139 


RLEO 


M8ine 


V 






13001+0527 


RTVIR 


MSm 


V 








14219+2555 


RXBOO 


M7.5 


V 




V 




13114-0232 


SW VIR 


M7111 


V 








L^UUj- /Odd 


(7 /\r o 


IVlD.Dili. 


V 




.1 
V 




20248-2825 


TMIC 


M7111 


V 








04566+5606 


TXCAM 


M8.5 


V 


V 






20038-2722 


V1943 SGR 


M7ni 


V 








04020-1551 


VERI 


M5/M6IV 


V 








18050-2213 


VXSGR 


M5/M6ni: 


V 








07209-2540 


VYCMA 


M3/M4II: 


V 


V 




V a 


13462-2807 


WHYA 


M7e 


V 


V 


V 


V 


01037+1219 


WXPSC/IRC +10 011 


M9: 


V 


V 






16011+4722 


XHER 


M8 


V 


V 






20075-6005 


XPAV 


M6/M7111:p 


V 






S-stars AGB 














19486+3247 


XCYG 


S 


V 


V 






04497+1410 


oORI 


S+WD 


V 








22196-4612 


;rGRU 


S5+G0V 


V 


V 


V 




09076+3110 


RS CNC 


M6111S 


V 




V 




20120-4433 


RZSGR 


S4,4 


V 








01159+7220 


S CAS 


Se 


V 








15492+4837 


ST HER 


M6s 


V 








19126-0708 


W AQL 


S6.6+FV 


V 


V 






07245+4605 


YLYN 


M5Ib-IIS 


V 
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Table 1. Continued 



IRAS name 


Identifier 


type 


PACS 


PACS 


SPIRE 


SPIRE Remarks 








Imag. 


Spec. 


Imag. 


Spec. 


C-rich AGB 


01144+6658 


AFGL 190 


C 


V 








18240+2326 


AFGL 2155 


C 


V 








23166+1655 


AFGL 3068 


c 


V 


V 


V 


V 


23320+4316 


AFGL 31 16 /IRC +40 540 


c 


V 


V 






00248+3518 


AQ AND 


Nvar 


V 




V 




19314-1629 


AQSGR 


CII 


V 








10131+3049 


C1T6 


c 


V 


V 


V 


V 


09452+1330 


CW LEO /IRC +10 216 


c 


V 


V 


V 


V a,b 


11331-1418 


HD 100764 


C-R. 


V 






15194-5115 


IRAS 15194-5115 


C 


V 


V 






02270-2619 


RFOR 


C 


V 








04573-1452 


RLEP 


Clle 


V 








01246-3248 


RSCL 


CII 


V 




V 


V 


20141-2128 


RTCAP 


CII 


V 








12544+6615 


RYDRA 


C 


V 








21358+7823 


S CEP 


CII 


V 




V 


c 


18476-0758 


SSCT 


CII 


V 








04459+6804 


ST CAM 


N 


V 








18306+3657 


TLYR 


C 


V 








19390+3229 


TTCYG 


N 


V 




V 


c 


03112-5730 


TWHOR 


CII 


V 








23438+0312 


TX PSC 


CII 


V 




V 




10329-3918 


U ANT 


N:var 


V 




V 




03374+6229 


UCAM 


N 


V 




V 




10350-1307 


UHYA 


CII 


V 








06331+3829 


UU AUR 


CII 


V 








19233+7627 


UXDRA 


CII 


J 

V 




J 

V 




15477+3943 


VCRB 


N 


V 






20396+4757 


VCYG 


N 


V 


V 






10491-2059 


VHYA 


C 


V 


V 


V 




17389-5742 


VPAV 


c+ 


V 








10416+6740 


VY UMA 


CII 


V 








05028+0106 


W ORI 


CII 


V 








05418-4628 


W PIC 


Nvar 


/ 

V 








15094-6953 


XTRA 


C 


V 








12427+4542 


YCVN 


Clab: 


V 




V 




21197-6956 


YPAV 


CII 


V 








Post-RSG 














19114+0002 


HD 179821 / AFGL 2343 


G5Ia 


V 


V 




V 


19244+1115 


IRC +10 420 


F8Ia 


V 


V 






RCrB 














15465+2818 


RCRB 


GOIab:pe 






V 
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Table 1. Continued 



IRAS name 


Identifier 


type 


PACS 


PACS 


SPIRE 


SPIRE 


Remarks 








Imag. 


Spec. 


Imag. 


Spec. 




Post^AGB 




AFGL 2688 


C-PAGB 


V 


V 


V 


V 




10215-5916 


AFGL 4106 


PAGE 








V 




04395+3601 


AFGL 618 


C-PAGB 


V 


V 


V 


V 




17150-3224 


AFGL 6815 


0-PAGB 


V 


V 








09371-1-1212 


Frosty Leo Neb. 


PAGB 




V 




V 




11385-5517 


HD 101584 


PAGE 


V 










17436-H5003 


HD 161796 


0-PAGB 


V 


V 


V 


V 




22272-1-5435 


HD 235858 


C-PAGB 


V 


V 








06176-1036 


HD 44179 / Red Rectangle 


C-PAGB 


V 


V 


V 


V 


b 


07134-1-1005 


HD 56126 


C-PAGB 


V 


V 


V 


V 




13428-6232 


IRAS 13428-6232 


PAGE 


V 


V 








16342-3814 


IRAS 16342-3814 


0-PAGB 


V 


V 








16594-4656 


IRAS 16594-4656 


C-PAGB 


V 










19343-1-2926 


M 1-92 


PAGE 


V 










07399-1435 


OH 231.8-1-4.2 


0-PAGB 


V 


V 




V 





PNe 


19327-1-3024 


ED -1-30 3639 


PN 




V 




17047-5650 


CPD -56 8032 


PN 






17347-3139 


GLMP 591 


PN 




V 




14562-5406 


Hen 2-113 


PN/WC 




V 




17423-1755 


Hen 3-1475 


PAGB/Be 




V 




10178-5958 


Hen 3-401 


PN/Be 




V 




09425-6040 


IRAS 09425-6040 


PAGE 




V 




16279-4757 


IRAS 16279-4757 


PAGB 




V 




220364-5306 


IRAS 22036-1-5306 


PN 


V 


V 




11119-1-5517 


NGC 3587 


C-PN 


V 


V 


17103-3702 


NGC 6302 


O-PNe 




V 




01391-1-5119 


NGC 650 


PN 


V 


V 


18021-1950 


NGC 6537 


O-PNe 




V 




175844-6638 


NGC 6543 


PN 


V 




V 


18517-1-3257 


NGC 6720 / Ring Nebula 


O-PNe 


V 




V 


19574-1-2234 


NGC 6853 


PN 


V 




V 




NGC 7027 


C-PNe 


V 


V 




22267-2102 


NGC 7293 / Helix Nebula 


PN 


V 




V 


21282-1-5050 


PNG093.9-00.1 


PN 




V 





Massive stars 



10541-6011 


AG Car 


LEV 


V 




G79.29-t-0.46 


LBV? 


V 


18184-1623 


HD 168625 


LBV? 


V 


10520-6010 


He 3-519 


Of/WN 


V 


10211-5922 


HRCar 


LEV 


V 


19092-1-1646 


M 1-67 


WN8 


V 




NGC 6888 


WN6 


V 


11065-6026 


Wra751 


LEV 


V 



SNe remnants 



Cas A/SNR 111.7-02.1 


V 


V 


V 


V d 


05314-^2200 Crab / SNR 184.6-05.8 


V 


V 


V 


V 


17276-2126 Kepler's SN / SNR 004.5-H06.8 


V 


V 


V 




SN 1181 / SNR 130.7-H03.1 


V 


V 






Tycho's SN / SNR 120.1-1-01.4 


V 


V 


V 


V 



Remarks: (a) PACS-i-SPIRE spectroscopy obtained as part of PV. (b) PACS-i-SPIRE imaging obtained as part of SDP (data are public), (c) PACS 
imaging obtained as part of SDP (data are public), (d) PACS -I-SPIRE imaging obtained as part of SDP. (e) Data taken in PACS-SPIRE parallel mode. 
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Table 2. Stellar fluxes (in Jy) in the apertures listed between parenthesis (in arcseconds) of the AGB and post-AGB stars observed 
per October 17th. Add 1342100000 to the AOR number quoted to obtain the Observation ID. Date refers to the observing date 
(format yy-mm-dd), when two are listed the first refers to the PACS observations, the second to the date of the SPIRE observations. 



Object 


PACS blue 


red 


SPIRE PSW 


PMW 


PLW 


AOR 


Date 


Remarks 




70 ^lm 


160 ^lm 


250 jum 


350 //m 


500 jum 








AFGL 190 


42.9 (33) 


6.26 (31) 








89183, 89184 


10-01-12 




AFGL 2155 


54.1 (35) 


8.57 (34) 








92777, 92778 


10-03-26 




AFGL 2688 


2047(110) 


488 (118) 


119(106) 


48.9 (140) 


18.2 (92) 


95837, 95838, 88168 


10-05-05, 09-12-16 




AFGL 278 


22.7 (29) 


3.30 (28) 








89840, 89841 


10-01-28 




AFGL 3068 


175 (40) 


27.6 (34) 


6.93 (56) 


3.32 (136) 


1.26 (92) 


88178, 88378, 88379 


09-12-21, 09-12-16 




AFGL 3116 


73.0 (46) 


12.6 (50) 








88490, 88491 


09-12-23 




AFGL 618 


953 (104) 


196 (86) 


46.1 (101) 


19.0 (87) 


8.91 (92) 


93132, 93133,91183 


10-03-31, 10-02-25 




ff ORI 


294 (130) 


35.5 (82) 


6.18(57) 


2.84 (68) 


1.56 (92) 


104435, 104436, 92099 


10-09-13, 10-03-11 


c 


AQ AND 


0.914 (14) 


0.391 (28) 








88488, 88489 


09-12-23 


a 


AQSGR 


2.05 (14) 


0.354 (28) 








95702, 95703 


10-04-30 




C1T6 






8.81 (72) 


3.77 (85) 


1.77 (92) 


106689 


10-10-17 




CWLEO 


3653 (244) 


610 (253) 


138 (211) 


57.7 (325) 


23.2 (276) 


86298, 86299, 86293 


09-10-25 




EP AQR 


26.6 (31) 


4.14 (31) 








95460, 95461 


10-04-23 


b 


HD 161796 


109 (42) 


14.5 (44) 


2.39 (45) 


0.731 (63) 


0.254 (92) 


88155,91127,91128 


10-02-24, 09-12-15 




HD 179821 


439 (63) 


62.0 (55) 








93490, 93491 


10-04-03 




HD 235858 


60.2 (42) 


7.98 (28) 








96757, 96758 


10-05-20 


c 


HD 44179 


139 (49) 


32.5 (55) 


9.51 (50) 


4.46 (67) 


2.46 (92) 


85549, 85550, 83682 


09-10-09, 09-09-12 




HD 56126 






1.21 (45) 


0.190 (92) 


0.388 (63) 


93007 


10-03-28 




IRAS 11385-5517 


207 (55) 


62.8 (62) 








102339, 102340 


10-08-09 




IRAS 15194-5115 


106 (48) 


22.9 (72) 








90247, 90248 


10-02-05 




IRAS 16594-4656 


82.0 (39) 


12.3 (28) 








93052, 93053 


10-03-30 




IRAS 19343-1-2926 


116(44) 


33.9 (46) 








91951, 91952 


10-03-10 




IRC -10 529 


77.9 (40) 


12.1 (44) 








96034, 96035 


10-05-09 




XCYG 


53.5 (48) 


8.49 (38) 








88320, 88321 


09-12-20 




//CEP 


100 (86) 


8.47 (39) 








91945, 91946 


10-03-10 


b 


NML CYG 


786 (82) 


109 (35) 








95485, 95486 


10-04-24 


c 


NML TAU 


202 (69) 


27.6 (67) 


5.99 (63) 


2.71 (119) 


0.876 (92) 


90343, 90344, 91180 


10-02-10, 10-02-25 


c 


OH 127.8-hO.O 


87.9 (41) 


9.90 (31) 








89181, 89182 


10-01-12 




OH 26.5-1-0.6 


317 (44) 


35.4 (28) 








91817, 91818 


10-03-09 


c 


«CET 


245 (99) 


35.2 (90) 


6.87 (63) 


2.71 (73) 


1.27 (92) 


90335, 90336, 89423 


10-02-09, 10-01-17 


c 


«ORl 


1.79(14) 


0.368 (28) 








91108,91109 


10-02-23 




;rGRU 


48.6 (56) 


6.50 (47) 


1.54 (45) 


0.687 (63) 


0.313 (92) 


96799, 96800,93791 


10-05-21, 10-04-05 


c 


RAQR 






1.03 (45) 


0.423 (63) 


0.226 (92) 


88166 


09-12-16 




RCAS 






2.90 (45) 


1.51 (63) 


0.917 (92) 


88578 


09-12-25 


c 


RDOR 


193 (79) 


32.7 (57) 


9.63 (62) 


4.73 (154) 


2.19(92) 


97685, 97686, 88164 


10-06-03, 09-12-16 




RFOR 


9.00 (21) 


1.38 (28) 








88502, 88503 


09-12-24 




RHYA 


49.5 (37) 


7.95 (32) 


2.88 (47) 


1.50 (85) 


0.841 (92) 


102231, 102232, 88662 


10-08-06, 09-12-27 




RLE? 


16.5 (28) 


2.51 (28) 








90304, 90305 


10-02-07 




RSCL 






1.63 (45) 


0.692 (63) 


0.303 (92) 


88657 


09-12-26 




RTCAP 


2.34(16) 


0.363 (28) 








96036, 96037 


10-05-09 




RT VIR 


23.4 (39) 


3.80 (34) 








101181, 101182 


10-07-24 




RXBOO 






1.44 (25) 


0.798 (63) 


0.402 (92) 


88158, 


09-12-16 




RYDRA 


4.61 (19) 


0.810 (28) 








88338, 88339 


09-12-21 




RZSGR 


4.19 (22) 


2.88 (31) 








93532, 93533 


10-04-03 


b 
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Table 2. Continued 



Object 


PACS blue 


red 


SPIRE PSW 


PMW 


PLW 


AOR 


Date 


Remarks 




70jum 


160 


250 jum 


350 iim 


500 /zm 








S CAS 


19.6 (26) 


2.64 (28) 








105040, 105041 


10-09-22 




S CEP 


15.7 (27) 


3.03 (28) 


1.09 (45) 


0.759 (63) 


0.410 (92) 


88342, 88343, 103604 


09-12-21, 10-08-24 


a 


ST CAM 


3.04 (16) 


0.477 (28) 








106044, 106045 


10-10-07 




ST HER 


8.68 (26) 


1.23 (28) 








88324, 88325 


09-12-20 




f APS 


30.2 (37) 


4.38 (33) 


1.05 (45) 


0.476 (63) 


0.174 (92) 


93048, 93049, 103637 


10-03-30, 10-08-24 


a 


TLYR 


2.92 (15) 


0.673 (28) 








95694, 95695 


10-04-30 




TMIC 


17.2 (26) 


2.90 (28) 








93036, 93037 


10-03-30 


c 


TTCYG 


0.725 (14) 


0.139 (14) 








93500, 93501 


10-04-03 


a 


TWHOR 


3.55 (18) 


0.537 (28) 








88373, 88374 


09-12-21 




TXPSC 


8.34 (43) 


1.61 (35) 








88344, 88345 


09-12-21 


b 


U ANT 


6.35 (14) 


2.51 (28) 


0.491 (25) 


0.210 (25) 


0.074 (25) 


88474, 88475, 88161 


09-12-23, 09-12-16 


d 


UCAM 


8.65 (21) 


1.69 (28) 


0.891 (45) 


0.515 (63) 


0.239 (92) 


90307, 90308, 103608 


10-02-07, 10-08-24 




UU AUR 


8.52 (26) 


1.44 (28) 








106332, 106333 


10-10-12 




UX DRA 


2.08 (17) 


0.481 (28) 


0.156 (45) 


0.110 (63) 


0.099 (92) 


88340, 88341, 103603 


09-12-21, 10-08-24 


a 


V CRB 


3.01 (17) 


0.580 (28) 








88326, 88327 


09-12-20 




V CYG 


28.1 (37) 


4.97 (28) 








88462, 88463 


09-12-22 


b 


VERl 


11.3 (27) 


1.55 (28) 








89943, 89944 


10-01-31 




VHYA 






2.74 (50) 


1.31 (63) 


0.552 (92) 


88160 


09-12-16 




VPAV 


5.18 (19) 


0.845 (28) 








104278, 104279 


10-09-10 




VXSGR 


179 (62) 


22.7 (36) 








91799, 91800 


10-03-09 




VYCMA 


1096 (102) 


137 (28) 








94070, 94071 


10-04-07 




W AQL 


70.5 (58) 


11.3 (59) 








94084, 94085 


10-04-08 


b 


W HYA 






7.50 (49) 


3.39 (67) 


1.66 (92) 


89519 


10-01-19 


c 


WORI 


8.07 (22) 


1.41 (28) 








90965, 90966 


10-02-22 




WPIC 


2.44 (14) 


0.575 (28) 








90974, 90975 


10-02-22 




WXPSC 


105 (50) 


19.0 (68) 








88486, 88487 


09-12-23 




XHER 


25.6 (47) 


4.07 (40) 








88322, 88323 


09-12-20 


b 


XTRA 


7.90 (25) 


1.30 (28) 








93050, 93051 


10-03-30 




YCVN 


10.4 (27) 


1.89 (28) 


0.630 (45) 


0.468 (63) 


0.235 (92) 


88330, 88331, 88157 


09-12-20, 09-12-15 


b 


YPAV 


2.94 (16) 


0.468 (28) 








88157 


09-12-15 





Remarks: Sources are listed alphabetically. Internal errors in the fluxes are small and of order 0.01 Jy (PACS) and 0.02 Jy (SPIRE). The errors in 
the fluxes are dominated by the flux calibration accuracy, currently estimated to be 10% for the PACS 70 ;um band and better than 15% at 160 yum, 
and 15% for all SPIRE bands. 



(a) Aperture includes only the central star; (b) Aperture includes (part of) the extended emission; (c) Some contamination by extended emission 
in the aperture or sky annulus; (d) Only PACS blue and red and SPIRE PSW are essentially free from emission of the detached shell. 



